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We derive the first four terms of an expansion in m 2 H /m% of the total 
Higgs cross section through gluon fusion. At NLO we demonstrate the 
excellent convergence of this series to the known result keeping the exact 
top mass dependence. At NNLO there is no known exact result, and our 
work represents a thorough quantitative investigation of the effects of finite 
top mass at this order. We discuss the applicability of our approach, and 
present numerical results for the LHC and Tcvatron. 

PACS numbers: 12.38.Bx,14.80.Bn 

1. Introduction 

The forthcoming LHC experiments are expected to elucidate the mech- 
anism of electroweak symmetry breaking. The most popular theoretical 
models of this phenomenon invoke a new scalar field called the Higgs boson, 
and the LHC has been designed with this firmly in mind. To discover the 
Higgs it is important to thoroughly understand its production and decay 
modes, as well as any relevant backgrounds. In this talk we describe recent 
work concerning the production of a Standard Model Higgs boson through 
gluon fusion. This process proceeds via a massive top quark, and as one 
can see from Fig.[T]the leading order diagrams are already one loop. At this 
order an exact analytical expression for the cross section is known, but at 
higher orders in a s analytic formulae retaining the exact top mass depen- 
dence are only partially known [H El [3] , although a numerical code known 
as HIGLU [I] is available, which evaluates the NLO cross section exactly. 
Fortunately, one finds that the cross section is very well reproduced by 
working in an effective theory in which the top quark is integrated out, and 
then weighting with the leading order mass dependence. Comparisons with 
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Fig. 1. Leading order diagram for the process gg — > iJ. 



HIGLU show that this procedure works extremely well up to NLO [5J [61 [7] • 
This is usually taken as sufficient justification for also using the effective 
theory approach at higher orders; the NNLO contributions [HJ [U [10] are 
known only in the effective theory. 

To test the heavy-top limit approximation, we work in the full theory, 
including the top quark, and evaluate the cross section at NNLO as an 
asymptotic series in 1/Mf. The technology to perform this expansion is 
well known, and automatised in the q2e/exp framework [llj . We assume 
Mi heavier than all other scales in the problem, and express all Feynman 
integrals as convolutions of massive vacuum integrals with at most three 
loops, and massless 3/4/5-point functions through 2/1/0 loops, respectively. 
We note that the purely virtual contributions at NNLO have recently been 
calculated by two separate groups [T2J Q3] • 

We are able to directly evaluate the single real emission contribution in 
terms of hypergeometric functions, which we then expand in e with the Hyp- 
Exp p3] package. For the double real emission part we perform an expansion 
in powers of (1 — x) and then integrate term by term. To demonstrate the 
cancellation of infrared poles we must then of course similarly expand the 
single real part. After renormalising the coupling a s (for which we adopt 
the MS scheme), top mass Mt and gluon wave function (both in the on- 
shell scheme) we are left with only 1/e 2 and 1/e infrared poles. These are 
absorbed in the PDFs as per the usual mass factorisation procedure. 

The description provided here is necessarily brief - we refer the reader to 
Ref. [15] for more details, full analytic results and comprehensive references. 



The total cross section is a convolution of the partonic cross section with 
non-perturbative PDFs. Equivalently, we can write it as an integral over the 
luminosity, 



where x is related to the partonic centre of mass energy by x = M^/s. 
The small x region of this integral therefore corresponds to high energy. 



2. Small-x Behaviour 
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For the LHC the partonic cms energy can range up to 14 TeV, and so 
our assumption that the top mass is the heaviest scale in the process clearly 
breaks down here. Fortunately the behaviour of the cross section in the limit 
x — ► is known exactly, i.e. keeping the full top mass dependence [16]. We 
can therefore improve our results, which are valid for large x, by smoothly 
matching them onto the small-x results of Ref . [16] . 
The small-x results are, 



= 3(7 C (1) + O(x) 



-9a C {2) lnx + c + O(x), (2) 



for the NLO and NNLO cross sections respectively, where the coefficients 
C« and C( 2 ) are given in Ref. [16] in the form of a numerical table, out 
of which we construct simple interpolating functions. The constant c was 
not determined, and we set it to zero. Our strategy for matching these 
expressions with our results is to construct functions which have the correct 
behaviour in each of the limits x — > and x — ► 1, up to some order N in an 
expansion in powers of (1 — x). We write 
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where &gg' N denotes the expansion of the partonic cross section around 



x 



1 through 0((1 - x) n ) 



3. Numerical Results 

In order to strictly test the heavy-top limit, we apply a consistent 1/Mt 
expansion to the partonic cross section, without factoring the LO mass de- 
pendence into the higher order terms. At NLO we therefore define, 

a^°(Mn = a S ag 6 Pg 5(l - x) + ^*$(A#). (4) 

In Fig. [2] we compare the ^-channel NLO cross section, evaluated ac- 
cording to Eq. (jl]), including successive terms in the 1/Mt expansion, with 
the exact result from the numerical code HIGLU. We observe excellent con- 
vergence of the series at both the LHC and Tevatron. Unfortunately the 
low-x behaviour of the other subprocesses (qg, qq and also the NNLO chan- 
nels qq and qq') is not known. However, their numerical contribution at 
NLO is small, at the level of a few percent in the case of qg and at the 
permille level for qq. 
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Fig. 2. Ratio of the gg induced component of the NLO hadronic cross 

section as obtained from Eq. to the value obtained from HIGLU jlj, 

when keeping successively higher orders in l/M t (decreasing dash-length 

corresponds to increasing order) ; the dotted line is the result obtained from 

( l 

the pure soft expansion & gg 
of Eq.©. 



(i),JV 



through order 1/M t without the matching 
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Fig. 3. Ratio of the gg induced component of the NNLO hadronic cross sec- 
tion as obtained from Eq. ([3]) to the heavy-top result of Eq. {T]) (decreasing 
dash-length corresponds to increasing order in 1/M t ); the dotted line is the 
result obtained from the pure soft expansion a g 2 J' N through order 1/M t 6 
without the matching of Eq. ((3]). 



The success of our approach at NLO means we can confidently apply it at 
higher orders. In Fig.[3]we compare the NNLO total cross section in our 1/Mt 
expansion approach with that obtained in the heavy top effective theory. 
The LO mass dependence is factored in only up to NLO, so that we can 
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examine the truly NNLO mass effects. We observe very good convergence 
towards the heavy-top results, which assures us of the high quality of the 
latter. The two results differ by less that 0.5%, which clearly justifies the 
use of the heavy-top effective theory so far in the literature, and also in 
future studies. 
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